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ABSTRACT 

A degenerated three-dimensional finite element based on the total 
Lagrangian, Incremental, formulation of a three-dimensional layered 
anisotropic median Is developed, and its use In the geometrically 
nonlinear, static as well as dynamic, analysis of layered conposite 
plates and shells Is demonstrated via several example oroblems. For 
conparlson purposes a two-dimensional finite element based on the 
Sanders shell theory with the von Karman (nonlinear) strains Is also 
presented. The elements have the following features: 

• Geometrically l inear and nonlinear analysis 

• Static and transient analyses 

• Natural vibration (linear) analyse*, 

• Plates and shell elements 

• Arbitrary loading and boundary conditions 

• Arbitrary lamination scheme and lamina properties 

The element can be used, with minor changes. In any existing general 
purpose programs. 

The 3-0 dimensional degenerated element has cociputational 
simplicity over a fully th^et -dimensional element, and the element 
accounts for full geometric nonlienarities In contrast to the 2- 
dimenslonal elements based on the Sanders shell theory. As demonstrated 


via m^erlcal examples, the deflections obtained by the 2-0 shell 
element deviate from those obtained by the 3-0 element for deep 
shells. Further, the 3-0 element can be used to model general shells 
that are not necessarily doubly -curved. For example, the twisted plates 
cannot be modeled using the 2-0 shell element. Of course, the 3-0 
degenerated element is coi^jtational ly more demanding than the 2-0 shell 
theory element for a given problem. In simmar) , the present 3-0 element 
is an efficient element for the analysis of layered conposite plates and 
shells undergoing large displacements and transient motion. 
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1.1 Motivation 

Conposite materials and reinforceo plastics are Increasingly used 
In aircraft, space vehicles automobiles, and pressure vessels. With the 
1ncr^»ased use of f Iber-rel nforced conposltes as structural elements, 
studies involving the thermomechanical behavior of shell conponents made 
cf conposltes are receiving considerable attention. Functional 
requirements and economic considerations of design have forced designers 
to use accurate but economical methods of determining stresses, natural 
frequencies, buckling loads, etc. 

An accurate prediction of the behavior of shell structures requires 
a realistic modeling of the actual geometry, material behavior, and 
kinematic description of the conponents. The partial differential 
equations describing the large-deflection behavior of anisotropic 
conposite shells of arbitrary geometry are not amenable to classical 
analytical methods. Consequently, nunerical and approximate methods must 
be used to predict desired design quantities (such as stresses, 
frequencies, and buckling loads). In the last two and a half decides 
the finite element method has emerged as the most powerful 
analysis method of structural analysis. 

The majority of the research papers in the open literature on 
shells is concerned with bending, vibrjtion, and buckling of isotropic 
shells. As conposite materials are making thei r way into many 
engineering structures, analyses of shells made of such materials become 
inportant. Further, with the increased application of advanced fiber 
conposites in jet engine fan or conpressor blades and high performance 
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aircraft, studies Involving transient response of coitposite shell 
structures are needed to assess the capability of these materials under 
dynamic loads. At the time the project was undertaken, a f Inite-element 
analysis of the nonlinear transient response of laminated anisotropic 
shells was not available. 

FI nite-e'ement analyses of shell structures In the past have used 
one of the three types of elements: 1. A 2-d1mens1onal (2-0) element 

based on a two-dimensional shell theory: 2. A 3-0 element based on 
three-dimensional elasticity theory; and 3. A 3-0 degenerated element 
derived frjm the 3-0 elasticity theory. The 2-0 shell theory Is derived 
form the three dimensional continuum field equations via, for example, 
an analytical Integra* i-jn through the thickness Is employed to reduce 
the theory to a two-dimensional theory. In doing so the static and 
kinematic resultants are defined and used to derive the equations. In 
contrast to the 2-0 shell theory, In the 3-0 degenerated element, the 
shell geometry and displacement fields are discretized from the outset 
In the sense of finite elements, and the element contains full geometric 
nonlinearity. The unavailability of a convenient, general nonlinear 
shell theory makes the 2-0 shell element restrictive In Its use. The 
nonlinearity Included In the 2-0 shell element Is that due to the von 
Karman strains. In which the products of the derivatives of the transverse 
deflection are neglected. In contrast to the 2-0 shell 

theory, no specific shell theory Is enployed In the 3-0 degenerated 
element; Instead, the geometry and the displacement fields are directly 
discretized and interpolated as In the analysis of continuum problems. 

The 2-0 elements based on shell theory are the most economical, followed 
by the 3-0 degenerated element. While the 3-0 element is most accurate. 


at least In theory, the conputatlonal ccst prohibits Its use In the 
nonlinear trans- it analysis of shells. 

The present study Is motivated by the lack of a f 1n1te>e1ement 
analysis wf geometrically nonlinear transient response of laminated 
anisotropic shells. The present study Involved the development of a 3-0 
degenerated shell element for the analysis of a layered anisotropic 
medlun, accounting for full geometric nonlinearity and transverse normal 
and shear stresses. The following literature survey provides a 
background for the present work. 

1.2 A Review of the Literature 

There exist a er of theories for layered anisotropic shells. 
Many of these theories were developed originally for thin shells, and 
are based on the K1 rchhoff-Love kinematic hypothesis that plane sections 
normal to the undeformed midsurface remain plane and normal after 
deformation. Excellent surveys of various shell theories can be found 
In the works of Naghdl [1] and Bert [2j. Here we review the literature 
on conposite shel Is. 

The first analysis that Incorporated the bending-stretching 
coupling (due to unsymmetric lamination In conposites) Is due to 
Ambartsjnyan [3,4]. In his analyses Ambartsuny an assimed that the 
Individual orthotropic layers were oriented such that the principal axes 
of material symmetry coincided with the principal coordinates of the 
shell reference surface. Thus, Ambartsunyan 's work dealt with what is 
now known as laminated orthotropic shells rather than laminated 
anisotropic shells; in laminated anisotropic shells the individual 
layers are, in general, anisotropic and the principal axes of material 
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symmetry of the Individual layers do not necessarily coincide with the 
principal coordinates of the shell. 

In 1962 Ocng, Pister and Taylor [5] fornulated a theory of thin 
shells laminated of anisotropic material. The theory Is an exte(V;,u<.n of 
the theory developed ‘'y Stavsky [6] for laminated anisotropic plates to 
Donnell's shallow shell theory, Cheng and Ho [7] presented an analysis 
of laminated anisotropic cylindrical shells using Flugge's shell 
theory. A first approximation theory for the unsymmetrlc deformation of 
nonhomogeneous , anisotropic, elastic cylindrical shells was derived by 
Wldera and Chung [8] by means of the asymptotic Integration of the 
elasticity equations. For a homogeneous. Isotropic material, the theory 
reduces to the Donnell's equations. 

A1 1 of the theories discussed above are based on the K1 rchhof f-Love 
hypotheses In which the transverse shear deformation Is neglected. The 
Love's first approximation theories are expected to yield sufficiently 
accurate results when (1) the lateral dimension to thickness ratio Is 
large; (11) the dynamic excitations are within the low-frequency range; 
(111) the material anisotropy Is not severe. However, application of 
such theories to layered anisotropic conposlte shells could lead to as 
much as 30X or more errors In deflections, stresses, or frequencies. 

For example, the thick -walled conposlte cylinders used for aircraft 
landing gears require a laminated thick -snell theory for their analysis. 

The effect of transverse shear deformation and transverse Isotropy, 
as well as thermal expansion through the shell thickness were considered 
by Gulatl and Essenberg [9] and Zukas and Vinson [10]. Gulatl and 
Essenberg [9] showed that the circumferential displacement conponents 
and the twist couple would arise due to the anisotropy and transverse 
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shet',' deformation. Whitney and Sun [11] developed a shear deformable 
theory for laminated cylindrical shells that Includes both transverse 
shear deformation and transverse normal strain as well as expanslonal 
strains. Recently, Widera and Logan [12,13] presented refined theories 
for nonhomogrneous anisotropic cylindrical shells. 

As far as the finite-element analysis of shells Is concerned, 
layered conposite shells have not received nearly as much attention as 
ordinary shells. The works of Dong [14] on statically-loaded 
orthotropic shell of revolution, Dong and Selna [15] on free vibration 
of the «;ame, Wilson and Parsons [16] on static axl symmetric loading of 
arbitrarily thick orthotrop'.c shells of revolution, and Schrlt and 
Monforton [17] on laminated anisotrjpic cylindrical shells are the only 
ones that considered the finite element method before the 1970's (note 
that the latter reference Is the only one that considered laminated 
anisotropic shells). In the 1970' s there was an Increased Interest In 
the finite-element analysis of bending and vibration of laminated 
anisotropic shells. Apparently the first finite-element application In 
laminated anisotropic shells of arbitrary geometry Is due to Thottpson 
and Bert [18]t who treated free (I.e., natural) vibration of general 
laminated anisotropic thin shells. Other finite-element analyses of 
layered anisotropic conposite shells include the works of Panda and 
Natarajan [19], Shivakumar and Krishna Murty [20], Rao[21], Seide and 
Chang [22], and Venkatesh and Rao [23]. Recently, Hsu, et al. [24] and 
Reddy [25] presented finite-element analyses of laminated thick 
cylindrical shells and laminated thick doubly curved shells, 
respectively . 

I 

} 
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All -'f the literature cited above is limited to the small 
displacement theory of shells. In the analysis of thin iiexlble 
conposlte shells one should take large deflections Into account. 

Because of the high modulus and high strength properties that conposltes 
have, structural conponents could undergo large deflections before they 
become Inelastic. Therefore, an accurate prediction of dlsolacements 
and stresses Is possible only when one accounts for the geometric 
nonl inearl ty . 

Finite-element analyses of the 1 arge-dl spl ace*'Kint theory are based 
on the principle of virtual work or the associated principle of 
stationary potential enery . Horrlgmoe and Bergan [26] presented 
classical variational principles for nonlinear problems by considering 
Incremental deformations of a contlnuun. A survey of various principles 
In Incremental form 1n different reference configurations, such as the 
total Lagranglan and the updated Lagranglan fornulatlon. Is presented by 
Wunderl 1ch[27], In the total Lagranglan description , all static and 
kinematic variables are referred to the initial configuration. In the 
updated Lagranglan description all variables are referred to the current 
configuration. Stricklin et al. [28] presented a survey of various 
formulations and solution procedures for nonlinear static and dynamic 
structural analysis. The f orirul ations Included are the pseudo-force 
method, total Lagranglan method, the updated Lagranglan method, and the 
convected coordinate method. The solution methods Included are the 
solution by direct minimization of total potential, Newton -Raphson and 
modified Newton-Raphson, and the first- and second-order self correcting 


methods. 
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The only large-deflection analysob of laminated co'jposite shells 
that can be found In the literature are the static analysis of Noor and 
Hartley [29] and Chang and Sawamiphakdl [30]. Noor and Hartley employed 
the shallow shell theory with transverse shear strains and geometric 
nonllnearltles to develop triangular and quadrilateral finite 
elements. Chang and Sawamiphakdl presented a formulation of the 3-0 
degenerated element for gecinetrlcally nonlinear analysis of laminated 
conposite shells. The formulation Is based on the updated Lagranglan 
description and It does not Include any numerical results for laminated 
shePs. 

From the review of the literature It Is clear that the 3-0 
degenerated element has not been exploited fully for geometrically 
nonlinear analysis of lam'nated anisotropic shells. Further, the 
transient analysis has not been reported In the literature. GbldeJ by 
these observations the present work was undertaken In the fall of 1981. 

1 .3 The Present Study 

The present study was undertaken to develop a finite-element 
analysis capability for the static and dynamic analysis of geometrically 
nonlinear theory of layered anisotropic shells. The 3-0 degenerated 
element with total Lagranglan description Is used to analyze various 
shell problems. 

Following this introduction, a description of the 2-0 shell element 
Is presented in Chapter 2. In Chapter 3 a detailed discussion of the 3- 
0 degenerated element is given. Application and cotiparison of the two 
elements are Illustrated via a nunber of shell problems. 
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Chapter II 

A SHEAR DEFORMABLE SHELL ELEMENT 
2.1 Governing Equations 

The basic equations of tnree-dlmenslonal elasticity theory can be 
sinpllfled for thin flexible bodies. A set of simplifying assunptlons 
that provide a reasonable description of the behavior of thin elastic 
shells are as folloMs: 

1. the thickness of the shell Is small conpared to the other 
dimensions; 

2. the transverse normal stress Is negligible; 

3. normals to the reference surface of the shell before 
deformation remain straight but not necessarily normal after 
deformation; 

4. the thickness-to-radlus of the shell Is assuned to be small 
conpared to unity; and 

5. In the second order terms, the derivatives of membrane 
displacements are small conpared to the derivatives of the 
transverse displacement. 

The shell under consideration is conposed of a finite number of 
orthotropic layer?? of uniform thickness, as shown in Fig. 2.1. An 
orthogonal curvilinear coordinate system (Cj,C 2 *C) is chosen such that 
the and ^ 2 " curves are lines of principal curvature on the 
midsurface C * 0, and C-curves are straight lines perpendicular to the 
surface C “ 0. A line element of the shell is given by (see [31]) 

(ds)^ = [(1 + C/Ri)ajd^j]^ + [(1 + C/R2)®2‘^^2^^ ^ 
where a. and (i * 1,2) are the surface metrics and radii of 
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curvature, respectively. In general, and are functions 
of only . 

For constant values of a^, 02 * R^ . and R 2 , the equations of motion 
are given by (with Cq ■ y (l/R^ - I/R 2 ) and dx^ ■ 1 see [25]) 

3N,. 0, 5u, ll^^l 


&N2 Q2 


2 

a u. 


■5?- ("6 - ‘o^s* * sJT * IT ■ Pi ::r * C’’2 ”’3 <17 * 


at 


1 1 ®^^2 

1 2 at'^ 


aq. 

ax2 

A 

^2 

- q) • 


aM, 

— a . n. 

■ P., 

7 ^ ^ 

'">1 

ax2 


3 

at^ 

*6 

aM2 

A 




ax2 ‘ '"2 

" ^3 

* 


1 . 1 3.*" “3 


1 2 at 


^2 
a u. 


1 2 at‘ 

1 . 1 V-. 


+ [p, + p, (f- ♦ ^)] —r (2.2) 
^ 1 2 at^ 


where u^ are the displacements of the reference surface along the and 
C axes, are the rotations of the reference surface with respect to 
the 5^ -axis, q Is the distributed load, and are the stress and 
moment resultants, and is the shear force resultant: 


L 

(N. ,MJ * E / 0.(1, C)dC . i » 1,2.6 

'^k-l 


L « ^k 

' /, K 

'k-i 


o^ dC 


1 « 4,5, 


(2.3) 


Here P.j are the Inertias 


n 
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^ ^ ^ k-1 4. , 


and (1 • 4.5) are shear correction factors. 

The strain-displacement relations for a large rotation (small 
strains) theory of shells are given by (see [32]) 


Cl • €i * C<i 


C 2 - €2 + C <2 


C 4 • C 4 

^5 * ^5 


where 


au, u, , 5u, u, , du, u, - 0 U 2 2 , 

0 _ 1. . _3 4. i r/ i _lvZ . / L + 

dXi Ri 2 '■^bXi Rj' ^axj Ri^ ^bXi^ •' 

5u, u, , bu, u, 5 bu^ u, - ^“ 12 , 

4 • -sr * ST * i C<5IT * ST’ * ‘»ST * ST’ * 'se’ ’ 


bu 
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2 2 
bu^ 


bUo bu, U, bu, 

0 L + 2 /_! ^ U /__3 ^ 

bx, 'bx 0 


Rl''bX2 
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s;) ^ > i^) 


bu, u, 
0 _ . . 3 _2 

^4 " bx^ ■ R 2 


bu, u, 
0 . . . 3 __1 

H bXi ■ Rj 


* bXj Ri ^bXi Rj 
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^♦2 . I ,*^'^2 . '^ 2 . 

^2 ^ ^ 57 ^ * ft”’ 

a^2 8U2 5 u^ 

"'s " axj ^ axj ' ‘"0 ^axj ' 


( 2 . 6 ) 


Invoking the fourth assimption, the contribution from the underlined 
terms in Eq. (2.6) can be neglected. The transverse normal strain and 
normal stresses are neglected in the present theory. 

The shell constitutive equations are given by 


N 


. - A, ,c. + B, .ic . 
i ij j ij j 


B, •e'" + 0, .<• 


IJ J 


IJ J 


(i.j - 1.2.6) 


% 4^4 ^ 45*^5 


" ^ 45®4 ^ 55^5 


(2.7) 


Here A^j, B,j and O^j (i.j « 1.2.6) denote the extensional.* floxural- 
extensional coupling, and flexural stiffnesses: 


(A 





S / (l.C.C^)dC. 

k-1 C, 'J 


( 2 . 8 ) 


(Ic ) 

where L denotes the total number of layers, and are the plane- 
stress reduced stiffnesses of the k-th lamina, referred to the shell 
axes. 


Equations (2.1)-(2.8) conpletely describe the dynamic equilibriun 
of a layered, anisotropic shell. These equations can be solved in 
closed form for the small displacement theory' of simple-supported. 


♦Quantities A^^, A^^ and Agg are actually shear sti f fnessess, not 

extensional ones. 
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cross-ply plates under sinusoidal distribution of the transverse load 
(see [33,34]). In order to solve practically inportant problems that 
involve other loadings, boundary conditions, geometry, and lamination 
schemes, one nust consider approximate method of analysis. In the 
following section, an 1 soparameteric finite-element analysis of the 
shell equations ( 2 . 1 )-( 2 . 8 ) is presented. 


2.2 Fi nite-EI ement Model 

A typical finite element is a doubly-curved shell element whose 
projections on the ^5 an isoparametric quadrilateral 

(0) 

element. Over the typical shell element Q , the displacements (u^.U 2 , 
‘^ 3 *^ 1 ’ * 2 ^ are interpolated by expressions of the form, 

^ i 

j*l 
^ i 

• z ,5^) , i • 1,2 (2.9) 

1 j.l 1 J A 

where <(«j are the interpolation functions, and u;j and (^;| are the nodal 
values of u^ and respectively. For a linear isoparametric element 
with nine nodes (N=9), this interpolation results in a stiffness matrix 
of order 45 by 45. 

Substitution of Eq . (2.9) into the variational forinjlatlon of Eq . 
(2.2) yields an element equation of the form (see Reddy [34]) 

[K]{a} + [M]fi} - (F) (2.10) 

where {a} « {(u^}, {U 2 }, {U 3 }, U^}, {♦ 211 ^. element 

stiffness and mass matrices, respectively, and |F} is the force 
vector. In the interest of brevity, the coefficients of mass and 
stiffness matrices are included in Appendix I. 
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To conplete the approximation, we should approximate the time 
derivatives In Eq. (2.10). Here we use the Newmark direct Integration 
scheme (the constant>average-acce1erat1on method). Use of the Nevxnark 
method to Eq.(2.10) yields (see Reddy [34]) 

where 

[K] • [K] ♦ a^[M], (F) - ♦ M(aoU^n * ^ 

■ l/(3At^) , a^ ■ a^At , a2 ■ - 1 . (2.12) 

Once tne solution {a} Is known at t„+^ ■ (n+l)At, the first and second 
derivatives (velocity and accelerations) of (a) at t,,^^ can be conputed 
from 

(2-'3) 

where 83 ■ (1 - a)At, and a^ ■ oAt. 

The element equations ( 2 . 11 ) can be assembled, boundary conditions 
can be Inposed, and the resulting equations can be solved at each time 
step using the Information known from the preceding tline step 
solution. At time t ■ 0, the Initial values of {a} , (a} , and (a) 
(obtained by solving Eq . (2.10) at t * 0) (>re used to Initiate the time 
marching scheme. 

In the present study the nine-node rectangular Isoparametric 
element was employed. Analogous to the shear deformable theory of 
layered conposite plates [34], the present theory can be recognized as a 
shear deformable theory derived from the classical shell theory by 
treating the slope-deflection relations * 0 ) as constraints, and 
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Including the constraints Into the variational fornulatlon of the shell 
equations by the penalty function method (see Reddy elements 

derived using such theory are very stiff (so-called locking Is observed) 
for thin shells, but yield good results for moderately thick shells. To 
overcome the locking phenomenon, the reduced Integration technique (see 
Zienklewicz, Taylor and Too [35]) nust be enployed In the evaluation of 
the stiffness coefficients associated with the shear ener^ terms (I.e., 
penalty terms). More specifically, the 2x2 Gauss rule nust be used for 
shear terms (1 .e. , those Involving and Agg) and the standard 

3x3 cjuss rule nust be used for the bending terms when the nine node 
quadratic Isoparametric element. 
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Chapter 1 1 1 

DEGENERATED THREE OWENSIONAL FINITE ELEMEN T 

3.1 Introduction 

The primary objectl.-* of this chapter is to review the formilatlon 
of equations governing geometrically nonlinear raotlon of a continuous 
medltjn. Due to the nature of the present manuscript, only necessary 
equations are presented. For additional details the reader Is referred 
to [36-40]. 

We describe the motion of a continuous boc^y In a Cartesian 

coord1nat<b system. The sinultaneous position of all material points 

(I.e., the configuration) of the bo<^ at time t Is denoted by C^, 

and C and C^ ^ denote the conf 1 guratl'^n? at reference time t ■ t„ and 
0 t+At 0 

time t + At, respectively (see Fig. 3.1). In the updated Lagranglan 
descrip'.lon all kinetic and kinematic variables are referred to tne 
current configuration at each time and load step. In the total 
Lagranglan description all dependent variables are referred to the 
reference configuration. The updated Lagranglan Is more suitable for 
motions that Involve very large distortions of the body (e.g., high- 
velocity Inpact). The total Lagranglan Is more convenient for motions 
that Involve only moderately large deformations. In the present study 
^he total Lagranglan formulation Is adopted. 

3.2 Fornulatlon of the Incremental Equations of Motion 

Here we present a derivation of the equillbt iun equations at 
different time steps using the total Lagranglan approach. The 
coordinates of a typical point In C^ Is denoted by '"x • ,^X 2 ,'”x 2 ) . 

The displacement of a particle at time t is given by 
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°x or 


^ -°x 

*1 i 


(3.1) 


The Increment of displacement during time t to t ♦At Is defined by 


t*/l t 

ui - u, 


(3.2) 


The principle of virtual displacements can be employed to write the 
equilibrium equations at any fixed time t. The principle, applied to 
the 1 arge>d1 spl acements case, can be expressed mathematically as 
fol lows: 


/ u, 5u, dV^./ S,J 6( e,j)dV 


t+Atc ,/t+-At 


/ ^*^T. 6u. dA^ 6u, dV, 

A 1 1 0 \ 

0 0 


1 °“1 "’o 


(3.3) 


where siwmatlon on repeated Indices Is Implied; V^, A^, and denote, 
respectively, a volume element, area element, and density In the Initial 
configuration, S^j are the cotrponents of the second Plolo-KIrchhoff 
stress tensor, c^j the conpof.ents of the Grten-Lagrangl an strain tensor, 
T^ the conponents of boundary stresses, and are the conponents of the 
body force vector. The superposed dots on u^ denote differentiation 
with respect to time, and 6 denotes the variational symbol. In writing 
Eq. (3.3) It Is assuned that Is relatec! to the displacement 
conponents by the kinematic relations 


t + A i ,t + A ^ t+At ^ t+At 

ij 2 ^ “1 ,j j,i m,i 


t*-it 




(3.4) 


where u< j ■ bu./ox.. The strain conponents . can he expressed In 

' .j I j 1 J 

terms of current strain and Incremental strain conponents: 


r 
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■1j 


»'°o« ( JUAur ; 

"j <'“t.j * S.t * 'vi Vj’ 


* I '“l,J * “j.l * '“i».l "«i.J * “ni.l * i “".t “m,J 


^ ®1j * 


(3.5) 


where e^j and n^j denote the linear and nonlinear Incremental strains. 
The stress conponents can be deconposed Into two parts: 

'"'Sj ■ ‘'iJ * 'u ''•*> 

where S^j Is the Incremental stress tensor. The Incremental stress 
conponents S^j are related to the Incremental Green >Lagrang«: strain 
conponents, c^j * e^j by the generalized Hooke's law: 


^1j " ^Ijkt'kA’ 


(3.7) 


where are the conponents of the elasticity tensor. Using Eq. 

(3.4)-(3.7), one can be express Eq. (3.3) In the alternate form 


Po ^ ^ “l^'^l ‘^''o ^ ^1jkl^®kJl*^j ^kJl*®1j^‘^''o 
0 ^0 

* ®'tj ■'•'o Si ‘"ij “''o <3-8) 

0 0 

where S>l Is the virtual work due to external loads. 


3.3 Finite-Element Forinjiatlon 
3.3.1 Geometry of the Element 

Consider a solid three-dimensional element shown In Fig. 3.2. The 
coordinates of a typical point in the element can be written as 


■ ^ '-I’top • , 1 , ¥ <«t>boU»n 


(3.9) 
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' ••• f .: ■:•/ 

where n Is the nunber of nodes, a'** the finite-element 

Interpolation (or shape) functions, which. In the element take the value 
of unity at node 1 and zero at all other nodes, and ^2 are the 
normalized curvilinear coordinates In the middle plane of tne shell, 
and C Is a linear coordinate In the thickness direction and x| , x^, and 
Xj are the global coordinates at node 1 . Here ^ 

to vary between -1 and + 1 . Now let (see Fig. 3 . 2 ) 


'3k 


^’'k^op 


- (x ) 

' k ^bottom 


“1 

e 


1 - 1 /I ’ I 

3 = ^3/1^3! 


( 3 . 10 ) 


where v^j^ Is the k-th conponent of the vector Vj. Then Eq. ( 3 . 9 ) 
becomes 


where h. Is the thickness of the element at node j. For small 

J 

deformation, the displacement of every poinc In the element can be 
written as 


J ' 


( 3 . 12 ) 


1 i ^ i ^ i 

where 0^ and 02 are the rotations about (local) unit vectors e^ and 62# 

respectively, Ui» U2» and U2 are the displacement conponents 

corresponding to the global coordinate x^ , X2, directions 

respectively, and u| , u^ and u^ are the values of the displacements 


; 
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(referred to x) at node 1. In writing Eq. (3.12), we assuned that a 
line that Is straight and normal to the middle surface before 
deformation Is still straight bu*” not nocessarl ly 'normal' to the middle 
surface after deformation. The strain ener^ corresponding to stress 
perpendicular to the middle surface Is Ignored to Inprove numerical 
conditioning when the three-dimensional element Is employed. This 
constraint corresponds only to a part of the usual assunptlons of a two- 
dimensional shell theory. The relaxation of the requirement that 
straight lines perpendicular to the middle surface remain normal to the 
deformed middle surface permits the shell to experience shear 
deformation - an Inportant feature In thick shell situations. 

3.3.2 Displacement Field In the Element 

In the present study the current coordinates ^x^ are Interpolated 
by the expression 



n 


E 

j-1 







t:j ) 
®31 > 


and the displacement by 


(3.13) 


‘“i *j - Ml” 


“i-/., (“‘M, -M,)’ 


(3.141 

(3.15) 


Here ^u!j and u"^ denote, respectively, the displacement and Incremental 
displacement conponents In the x^-direct1on at the j-th node. The unit 

«i 

vectors e, and can be obtained from the relations 


r 


! 
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■ ‘ j 

rcY 



t:i 


53" 



(3.16) 


where is the unit vector along the (global) X2*ax1s. If we assume 
that the angles e| and are very small, then we can write 


*1 

53 


ri-i ^ t“i.i 

- i2»l * 5|«2 


(3.17) 


Substituting Eq. (3.17) Into Eq. (3.15), we obtain 


u^ « Z (PjCu-j + jChj(- + ^e^^0^)] 

J ^ 


(3.18a) 


or 

fu} « [T](a} (3.18b) 

where fu} Is the column of three displacements at a point, {a} Is the 

column of 5n (five per node) displacements: u^ , 0^, J * ^*2 '’J ^ 

■ 1,2,3, and [T] Is the transformation matrix defined by Eq. (3.18a). 
Thus, for each time step one can find the normal vectors from Eq.(3.16) 
and (3.17) and the Incremental displacements at each point from Eq. 
(3.18), once the five generalized displacements at each node are 
known. Next we discuss the procedure of determining the generalized 
displacements of an element. 

3.3.3 Element Stiffness Matrix 

The strain-displacement equations (3.4) can be expressed in the 
operator form 


f 


t 


r 
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(3.19) 


(e| • CA]{u^} 

where (ej ■ (e^^ «22 ®33 ^*12 ^*13 ^*23^^’ * function 

of qU^ j, and {u^} Is the vector of the conponents of the displacement 
gradient 



The vector (u^| Is related to the displacement Increments by 

[%] - CN]{u} • [N][T]{a} (3.20) 

and 

ie} - [A][N]CT](a} 

5 CB](a} (3.21) 

where [N] Is the operator of differentials. 

Substitution of Eq. (3.21) Into Eq. (3.8) yields 

/ (^D-]^u}dV^ ♦ + ^DC,<l]){a} • (3.22) 

0 

where (R}, and {F} are the linear and nonlinear stiffness 

matrices, force vector, and unbalanced force vectors: 

■ / ‘CB]^CC] ‘CBJJVj, . . / ‘[B]^CS] 

^0 ^0 

(F| • / ‘fS]^(SidV„ (3.23) 

V 

0 

A 

Here [S] and (S} denote the matrix and vector, respectively, of the 
second Plola-Ki rchhoff stress. 

Since we are dealing with laminated conposite structures, the 
Ifiportant thing Is how to perform the integration through the 
thickness. One way Is to pick Gaussian points through the thickness and 
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then no explicit Integration through the thickness Is performed. The 
CPU time will be Increased If the nimber of layers Is Increased, because 
the Integration should be performed separately for each layer. The 
other way Is to perform explicit Integration through the thickness and 
reduce the Integral to a two-dimensional problem. The Jacobian matrix. 
In general. Is a function of Cp Zienklewicz, et a1. [35] 

suggested that terms In c to the first power may be neglected, provided 
the thickness-to-curvature ratios are small. This approximation Implies 
that the derivatives of x^ with respect to Cp C 2 * ^ 

substantially the same at either end of a ml d- surface-normal line. Thus, 
the Jacobian [J] becomes Independent of C end explicit integration can 
be enployed. If C terms are retained In [J], Gaussian points through 
the thickness should be added. In the present study. It Is assigned that 
the Jacobian Is Independent of C* 

3.3.4 Time Integration and Mass Matrix 

Any attenpt to solve Eq. (3.22), whether by direct Integration or 
by modal analysis, should take advantage of the symmetry and bandedness 
of the stiffness matrix. The Initial conditions for Eq. (3.22) are the 
displacements and velocities at time t ■ 0; therefore, the direct 
Integration requires the Information at the previous time t. In order to 
predict the state of motion at the current time t+6t. The direct 
Integration techniques can be divided Into two types: explicit and 

InpHclt Integrations (see [41-45]). Inexplicit Integration, we 
solve u at time t+6t based on the equlllbriim conditions of the 
Structure at time t. The central difference method Is an example of 
explicit Integration method. In the implicit method the solution at 
time t+flt Is based on the equillbriun condition of the structure at 
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time The Houbolt, Wilson, and Nevmiark methods provide examples of 

the Implicit method. In the present study the Netenark direct 
Integration scheme Is enployed. 

The NeMiiark Integration scheme can be thought of as an extension of 
the linear acceleration method: 

• Mil * A^lil ♦ C(j - e)‘lM ‘ 

-'ui ‘CO -Y)‘iA} *,““u)]At 

where {a} Is the generalized displacement vector of any point 

and 6 and y are the dimensionless parameters of generalized 

acceleration. Chan et a1. [46] have discussed the special case 

^ " IT ^ “ T* coincides with a procedure developed by Fox and 

Goodvin [47]. For the constant average acceleration we have 

P ■ j and Y ■ y. and for the linear acceleration P ■ and Y ■ "j* 

To apply the Nevmiark Integration scheme to the equlllbrlim 
equations (3.22), we start from 

t+At|^}(k) ^ t+At|^|(k-1) ^ 25) 


where k Is the Iteration nunber. The velocity and acceleration of any 
point In the element can be written as 


t 


1 


n ^ . n . • . 

E 4<i uj + E T 4»iCh.(’^e^. 

j«l J ’ j«l ^ J J 



) 


(3.26) 


r 


••''iMTi p^r?r 13 
Cr . » (. ‘ 7Y 


” f*i ^ 1 f^i 

E 4». u , ■♦• E w 4>iCh i( 

j.l J J 1 «1 ^ J J 


o:j ) 
*3r 


(3.27) 


Since we are dealing with transient problerc, the load vector ( R} can be 
a function of time. From Eq . (3.27) we have for the 1-th conponent 


^1 *o' ^1 " ^1 ' *2 ^1 *3 ^1 


tj ^ . f; . , t+At" 
^1 “ ^1 *4 ^1 * *5 ^1 


(3.28) 


where 


° S(At)^ 


'1 pAt ‘ ^2 pAt 


*3 “ Ip " ^ . a^ « At(l - y) . «5 ■ yAt 


Substituting Eq. (3.22) Into Eq.(3.28), we obtain 


(3.29) 


• a (ejpj - ) - a - a ^e^ + ^e^ 

®31 *o^®l®21 ®2®ir ^2®31 ®3 ®1 ®31 


. t-j t;j mt;j 

®31 ®31 ®4 ®31 ®5 ®31 


(3.30) 


Finally the conplete approximation of the equations of motion over an 
element becomes 


(a^’^CM] + nK])fA^'^h - - t"-^t|p(k-l)| 
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- It - 'iPjD) * ajIPil (3.31) 

Mhere 

CH) ■ / p„ 'cT] dv„ 

y U U 

0 

(Pll ■ P„ '|A| [T]dv„ 

IPjl • / P„ **^'|A|"‘'”CT]dV^ 

0 

P„‘|Al[T]dV^ 

0 

IV ■ / Po '|a| CT]dV^ (3.32) 

^0 

Equation (3.31) can be expressed in the following final form: 

[K](A| • - t*At(p|(k-l) I3 J3J 

where 

[k] • A, ‘cm: ♦ ‘[K] 

..j c'lfil - it '''"al - 'I V' > 5 * ^3 1 " 4 I 

(3.34) 

This completes the finite-element formulation of the 3-0 degenerated 
element. 
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CHAPTER IV 

NIHERICAL VALIDATION OF THE ELEMENTS 

4.1 Introduction 

The present chapter is devoted to the validation of the finite 
elements developed herein based on the two-dimensional shell theory and 
three-dimensional contlnuim theory. The elements are validated by 
conparing the present results with those available In the literature for 
Static bending, natural vibration, and transient response of Isotropic 
plates and shells. Nunerical results are presented to bring out the 
limitations and restrictions of the present elements. All of the 
results presented here were obtained on an IBM 370/3081 conputer with 
double precision arithmetic. 

The results to be discussed are grouped Into three major 
categories: (1) static bending, (2) natural vibration, and (3) 
transient response. All results, except for the vibrations, are 
presented In a graphical form. 

4.2 Static Analysis 

Here we present a discussion of four exanple problems, all 
involving shel 1 structures. 

4.2.1 Cylindrical Panel Under the Influence of Gravity (I.e., under its 
own weight ) 

Consider the circular cylindrical panel shown in Fig. 4.1. The 
geometric parameters and material properties are listed below: 

R » 25.0 ft E » 3 X 10® psi 

a * 25.0 ft V » 0.0 

h » 3.0 in. g * 90 psi 


9 . 40* 
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supported 
by a rigid 

diaphram: u 


Figure 4.1 



Geometry of the cylindrical shell used In Problem 1 
of Section 4.1 
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The panel Is supported on rigid diaphragms on the curved edges and free 
on the straight edges. We wish to find the deformation of the panel 
under Its own weight. Many authors nave enployed this problem as a 
standard test problem for checking new elements or numerical schemes. 

The radial displacements along the central curved line are shown In Fig. 

4.2 for both the 3>0 element and the 2-0 deformable shell element. The 
longitudinal displacements along the supported edges are shown In Fig. 
4.3. A uniform mesh of 2 x 2 eight *node elements vas used In one 
quadrant for this analysis. The solutions agree ve y closely with those 
obtained in Reference [48]. 

4.2.2 Cylindrical Shell Subjected to Radial Pressure 

Consider the circular cylindrical panel shown 1n Fig. 4.4. The 
shell is clatiped along four edges and subjected to uniform radial Inward 
pressure. The loading Is nonconservative, that Is, the direction of the 
applied load Is normal to the cylindrical surface at any time during the 
deformation. The geometric and material properties are 
R > 2540 mm, a * b ■ 254 mm, h ■ 3.175 mm, 

9 - 0.1 rad, E ■ 3.10275 k.N/mm^, v ■ 0.3 
Due to the symmetry of the geometr 7 and deformation, only one quarter of 
the panel Is analyzed. A load step of 0.5 KN/m^ was used In order to 
get a close representation of the deformation path. Fig. 4.5 shows the 

t 

central deflection versus the pressure for the first of the three sets 
of panel dimensions: 

(1) a ■ 254 mm, b » 254 mm ( deep shell) 

(11) a » 635 mm, b * 635 mm 

(111) a ■ 1270 mm, b » 1270 mm ( shallow shell) 

In cases (11) and (111) the analysis was limited to linear solutions. 


in ft.) 
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0. 2. 4. 6. 8. 10. 12. 

Center deflection, w (In mm) 

; 

Figure 4.5 Load-deflection curve for the clamped cylindrical shell , 

problem discussed In Section 4.2.2 i 




36 


Table 4.1 contains a conparison of the center deflections obtained by 
using the two elements for various meshes. From the results presanted 
in the table it is clear that for shallow shells both elements give 
almost the same results. The difference in the deflections predicted by 
the two elements increases as the shell becomes shallow (case (iii)). 

4.2.3 Cylindrical Shell Subjected to Center Point Load 

Figure 4.6 shows a circular cylindrical shell which has the same 
geometric (except, h ■ 12.7 mm) and material properties ac the one in 
Problem 2. The longitudinal boundaries are h^;.g,.'d and immovable, and 
the curved edges are free. A concentrated pcint load is applied at the 
center. One quarter of the panel w ,s analyzed using a 2x2 mesh of nine- 
node quadri lateral elements. The load step is 0.5kN. Figure 4.7 
contains the plot of the central deflection versus the load. The 
results agree very closely with those obtained by Ohatt [49]. 

4.2.4 Spherical Shell Subjected to Point Load at the Center 
The spherical shell shown in Fig. 4.8 is subjected to a 

concentrated load P at the crown. The boundaries are all hinged and 
immovable. The geometric parameters and material properties are given 
below: 

* ^2 * 2540 mm 
a * L) * 784.9 mm 
h » 99.45 mm 
E ■ 68.95 N/mm 
V * P.3 

One quarter of the shell was analyzed by a 2x2 mesh of nine node 
quadrilateral elements. The load step is 10 kN. Figure 4.9 contains a 
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Table 4,1 Conparlson of linear center deflections obtained by the 2-0 
and 3-0 elements for Problem 2. 


Case 

Mesh 

3-0 

2-0 


2x2 

0.52579 

0.52321 

(1) 


3x3 

0.52505 

0.52265 


2x2 

0.35888 

0.39431 

(11) 


3x3 

0.3205 

0.34228 


2x2 

0.30636 

0.37667 


(ill) 


3x3 


0.28674 


0.27021 
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Figure 4.6 Geometry of the cylindrical shell problem 
discussed in Section 4.2.3 




P 


Center transverse deflection, w (In an) 


0' ^CLi,' 



Figure 4.7 Load-def1ect1on curve for the cylindrical shell 
problem discussed in Section 4.2.3 


all edges are hinged; 
u « w • 1 /,^ « 0 on BC 

* i'x “ ^ on AB 


“ »<2 “ 2540 (im, 
b » 784.9 mn, 
99.45 mm. 


Figure 4.8 Geometry of the spherical shell discussed 

Section 4.2.4. 
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conparlson of the center deflections obtained using the 2-0 and 3-0 
elements, which agree with that of Leicester [50], who used a shallow 
shell 2-0 finite element. 


4.3 Natural Vibration of Cantilevered Twisted Plates 

Here we discuss the results obtained for natural frequencies of 
cantilevered twisted plates. This analysis was motivated by their 
relevance to natural vibrations of turbine blades. Consider a 
cantilevered plate with a twist angle 9 at the free end. The plate is 
made of an isotropic material. Table 4.2 contains the natural 
frequencies of a square plate for various values of the twist angle 9. 

A 2x2 mesh and 4x4 mesh of nine-node elements are employed to study the 
convergence trend. The results of the refined mesh are included in the 
parentheses. The results agree with many others published in a recent 
NASA report. Tables 4. 2-4. 5 contain natural frequencies of twisted 
plates for various aspect ratios and si de-to-t hick ness ratios. 

4.4 Transient Analysis 

Here we present results of the nonlinear transient analysis of a 
cantilevered beam and a spherical shell. Both problems have been solved 
by other investigators. 


4.4.1 Cantilevered Beam Under Uniformly Distributed Load 

Consider a cantilever beam under a uniformly distributed transverse 
step load, as shown in Fig. 4.10. The geometric parameters and material 
properties are given below: 


OR'Cr 

Of f c-:;'- 



L * 1 0 i n 
h » 1 1 n 
b » 1 i n 


1 

) 

i 


r 


Center deflection, w (in mn) 
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Table 4.2 Natural Frequencies of Twisted Plate Vibration (a/b • 1, a/h • 20) 


w ■ ua^pl\/0 , 0 ■ ■ ■ , V ■ 0,3 

12(lV) 


Twist 

Angle 

1 

2 

Hade 

3 4 

5 

6 

0* ** 

* 3.4556 

**(3.4583) 

8.4110 

(8.3353) 

22.0999 

(21.0238) 

28.2089 

(26.7465) 

31 .9740 
(30.1454) 

55.1625 

(52.0784) 

15* 

3.4359 

10.2920 

21.5199 

27.2054 

32.7430 

44.5375 

30* 

3.3790 
(3 J694) 

13.7014 

(14.2222) 

19.9840 

(18.9795) 

25.0943 

(26.8104) 

34.3341 

(34.4591) 

45.8987 

(45.7547) 

45* 

3.2908 

18.1009 

15.9097 

23.5680 

35.5332 

45.7013 

60* 

6.1800 

17.8319 

15.5635 

24.1842 

36.1466 

44.9152 


*2x2 ,9 -node mesh 

**4 X 4,9-node mesh 
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Table 4.3 Natural Frequencies of Twisted Plate Vibration (b/a ■ 3, a/h ■ 20, 
3x3, 9-node mesh) 


u ■ ub^/ph/D , D ■ — — 5— , V ■ 0.3 

12(l-v^) 


Twist 

Angle 

1 

2 

3 

Mode 

4 

5 

6 

7 

0* 

3.4150 

20.8772 

21.6190 

65.9706 

66.2590 


127.256 

15* 

3.4009 

20.8798 

22.1 118 

21.6032 

68.0938 

69.3258 

130.284 

30* 

3.3598 

19.4048 

25.3743 

60.2183 

73.5180 

7 7.4493 

138.176 

cn 

0 

3.2956 

17.5289 

29.8404 

58.2600 

80.9488 

88.5246 

148.8975 

60* 

3.2136 

15.7431 

34.8827 

55.8921 

89.2028 

100.7760 

155.070 
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Table 

4.4 Natural Frequenc1e< 
u ■ ub^/ph/0 , 0 

of Twisted 

Plate Vibration (a/b 
, V ■ 0.3 

■ 1, a/h ■ 5) 

12(l-v^) 

Twist 



K)de 



Angle 

1 

2 

3 

4 

5 

6 

0* ** 

* 3.33916 

7 .3948 

10.8083 

18.4930 

23.7907 

26.0552 


**(3.3390) 

(7.3559) 

(10.883) 

(17.757) 

(22.769) 

(24.125) 

15* 

3.31713 

7.4816 

10.8053 

18.4043 

23.6767 

24.9474 


(3.3170) 

(7.4504) 

(10.774) 

(17.771) 

(22.694) 

(24.083) 

30* 

3.2538 

7.7593 

10.5248 

18.4091 

23.3734 

24.61 16 


(3.2538) 

(7.7089) 

(10.478) 

(17.795) 

(22.471 ) 

(23.943) 

45* 

3.1570 

8.1435 

10.1270 

18.3843 

22.9126 

24.0566 


(3.1569) 

(8.0728) 

(10.062) 

(17.79) 

(22.1 17) 

(23.651 ) 

60* 

3.0370 

8.5855 

9.67198 

18.3089 

22.3670 

23.3533 


(3.0366) 

(8.4814) 

(8.5911) 

(17.730) 

(21.684) 

(23.160) 


* 2 X 2.9-node mesh 

**3 X 3,9-node mesh 
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Table 4.5 Hatural Frequencies of Twisted Plate Vibration (b/a ■ 3, a/h ■ 5, 
3x3 ,9-noce nesh ) 



tal ■ ub^ 

/ph/b , 0 ■ 

Eh-* 

■ 0 ^ 




W— , V 

12(lV) 

■ U • J 



Twist 

angle 

1 

2 

Mode 

3 

4 

5 

6 

0* 

3.3908 

15.551 

19.124 

21.065 

59.924 

61 .949 

15* 

3.31 61 

15.192 

19.231 

21.572 

60 .088 

60.830 

30* 

3.3336 

14.379 

19.549 

22.81 1 

60.576 

58.472 

45* 

3.2674 

13.449 

20.060 

24.404 

61 .360 

55.874 

60“ 

3.1833 

12.548 

20.741 

26.139 

62.416 

53.381 
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E ■ 1.2 X 10® psi 

V • 0.2 

p ■ 10 lb sec^/i 

The cantilever Is analyzed using 'our nine-node 2-0 shell and 3 0 
degenerated elements (see Fig. 4.11). The load applied Is 
nonconservative, because the load follovs the deformed beam and stays 
normal to It at all times. Figure 4.11 shows the tip deflection versus 
load obtained In the static analysis. Figure 4.12 contains plots of tip 
deflection versus time; obtained using the 2-0 and 3-0 elements and by 
Bathe et a1. [37]. The time step enployed Is 1.35 x 10"^ sec. The 
reason for the difference between the solutions predicted by the shell 
element and 3-0 element Is that the load In the 3-0 element follows the 
deformed shape and Is perpendicular to the deformed beam, whereas In the 
2-0 shell theory It Is always vertical. Consequently, the vertical load 
conponent Is larger In the 2-0 shell element than In the 3-0 element, 
and thus explains the difference In the solutions. 

4.4.2 Spherical Cap Under f lsymmetrlc Pressure Loading 

Consider a spherical cap, clanped on the boundary and subjected to 
axisymmetric pressure loading. The geometric and material properties 
are 

R * 22.27 in 
h « 0.41 m 
E « 10.5 X 10® psi 
V « 0.3 

pg - 0.095 1b/1n^ 

0 » 26.670 
p^ » 100 psi 
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^ ■ lO"^ sec 

This problem has been analyzed by Stricklin, et a1. [51] using an 
axisymmetrl c shell element. In the present study the spherical cap Is 
discretized Into five nine-node 2-0 and 3-0 elements. Fig. 4.13 shows 
the center deflection versus time. The present solutions are In 
excellent agreement In most places with that of Stricklin et a1. [51]. 
The difference between the solutions Is mostly In the regions of local 
minima and maxima. 


r 


I 

I 
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Figure 4.13 Center transverse Jisplacement versus time 

for a spherical cap under axi symmetri c dynamic 
loading ( load « 100 psi.) 
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CHAPTER V 

NUMERICAL RESULTS FOR CO»f>OSITE PLATES AND SHELLS 

5.1 Introduction 

In this chapter we discuss the numerical results obtained by the 2* 
0 and 3-0 elements for the nonlinear analysis of layered anisotropic 
conposlte plates and shells. Most of the results presented here are na< 
and therefore cannot be conpared with other results to make quantitative 
Judgements concerning their accuracy. Results for both static bending 
and transient •'esponse are discussed. 

5.2 Static Analysis 

5.2.1 Orthotropic Cylinder Subjected to Internal Pressure 

Consider a clamped orthotropic (E2 ■ 2.0 x 10® psi, E^/^2 “ 3.75, 
612/E2 ■ 0.625, V ■ 0.25) cylinder of radius R ■ 20 In and length 20 In, 
and subjected to Internal pressure Pq ■ 6.41 n psi (see Fig. 5.1). The 
problem was analyzed, for linear deflections only, by Rao [21] using 
shallow thin shell elements. A mesh of 2x2 nine-node elements Is used 
In an attempt to analyze the problem. The linear center deflections 
obtained by the 2-0 and 3-D elements are 0.0003764 In., and 0.0003739 
In., respectively. These values coirpare favorably with 0.000366 In. of 
Rao [21] and 0.000367 of Timoshenko's analytical solution [52]. The 
latter two solutions are based on classical shell theory. 

In the large-deflection analysis the present results are conpared 
with those of Chang and Sawamiphakdl [30]. A value of 2.5 ksl Is used 
for the load step. Figure 5.2 shows a conparison of the present 
deflection with that of [30], which used a 3-0 degenerated element based 
on an updated Lagrangian approach. The agreement is very good. 
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Internal pressure (ksi) 
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Figure 5.2 Center transverse deflection versus 
internal pressure 


p 
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5.2.2 N1 ne-Layer Cross-Ply Spherical Shell Subjected to Uniform Loading 
Consider a spherical shell cross-ply laminated of nine layers of 

graph1te-epo)yr material (E 1 /E 2 * ®23^^2 “ ®13 “ ^12 " ^2 ’ 

and v ^2 * 0*3) » subjected to uniformly distributed loading, and sinply 
supported on all Its edges (I.e., transverse deflection and tangential 
rotations are zero). In Fig. 5.3, a conparlson of the load*<Jef1ect1on 
curves obtained by the present elements with those obtained by Noor and 
Hartley [29] Is presented (for the parameters h/a ■ 0.01 and R/a ■ 

10). The results agree very well with each other, the present 2-0 
results being closer to Noor and Hartley's solution. This Is expected 
because their element Is based on a shell theory. 

5.2.3 Two-Layer Cross-Ply and Angle-Ply (45* /-45*) Cylindrical Shells 

Uhdrr'thifoVm' Loading — 

The geometry of the circular cylindrical shell used here Is the 

same as that shown In Fig. 4.1. The shell Is assumed to be simply 

supported on all It; edges. The material properties of Individual 

lamina are the same as those used In Problem 2 of this chapter. A mesh 

of 2x2 nine-node elements In a quarter shell Is used to model the 

problem. The results of the analysis are presented In the form of load- 

deflection curves In Fig. 5.4. Frcm the results, one can conclude that 

the angle-ply shell is stiffer than the cross-ply shell. This can be due 
to the bending-stretching coupling. 

5.2.4 Two-Layer Cross-Ply and Angle-Ply (45**/-45°) Spherical Shells 
Under Uniform Loading 

The geometry and boundary conditions used in this problem are the 
same as those used In Problem 2 of this chapter. The geometric 
parameters used are: R/a » 10, a/h » 100. The load-deflection curves 

for the cross-ply and angle-ply shells are shown In Fig. 5.5. From the 


Center deflection, w/h 
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Figure 5.3 Deflection versus load parameter for nine- 
layer cross-ply (0V90V0*/. . . ) spherical 
shell discussed in Section 5.2.2. 


1 

1 

I 

I 


Nondl mens Iona 11 zed center deflectlont w/h 


r- 



Figure 5.4 Deflection versus the load parameter for two-layer 
composite cylindrical shell (see Figure 4.1) 
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AP 

0.24 

0.20 

0 . 16 , 

0 . 12 , 

0 . 08 , 
0.04 • 



Figure 5.5 Nondimenslonal ized deflection versus; tho 
for l«.(„ated shell, dl'L„ed 
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plot It Is apparent that, for the load range considered, the angle-ply 

shell, being stiffen, does not exhibit nuch geometric nonlinearity. The 

load-deflection curve of the cross-ply shell exhibits a varying degree 

of nonlinearity with the load. For load values between 100 and 150, the 

shell becomes relatively more flexible, (which can be due to bending- stret- 
ching coupling) 

5.3 Transient Analysis 

5.3.1 Two-Layer Cross-Ply Plate Under Uni form Load 

The problem Is the same as that analyzed by Red<1y [53] using a 
plate element. The geometric and material parameters used are 
a ■ b ■ 244 In., h ■ 0.635 In., ■ 0.5 x 10"^ psi 
Ej ■ 17.578 X 10® psi, ■ 0.7031 x 10® psi, ■ 0.5 x 10® psi 
p ■ 0.2547 X 10’® 1b-sec^/1n^, v ^2 ■ 0.25, 6t ■ 0.002 sec. 

Figure 5.6 contains a plot of the center deflection versus time obtained 
by the 3-0 element. The solution Is In excellent agreement with that of 
Red(^ [53]. 


5.3.2 Two-Layer Cross-Ply Cylindrical Shell Under Uniform Load 

A cylindrical shell with a ■ b ■ 5 In, R ■ 10 In, h ■ 0.1 In Is 
simply -supported on the four edges. The deep shell Is laminated with 


two layers (0*/90*) and loaded by a uniform step load P 


A 

E,h^ 


50. 


Figure 5.7 contains a plot of the center deflection versus time for 2-0 

-4 

and 3-0 elements with 6t ■ 0.1 x 10 sec. The solutions obtained using 
the two elements are In good agreement. In Fig. 5.8, the solid line 

A 

Indicates the center deflection versus time for load P ■ 1000 and time 

-5 

step at • 0.3 X 10 sec. The amplitude Is almost twelve times that due 

A 

to load P ■ 50, v/hereas the load Increases twenty times. The dotted 


t 
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Center deflection, w (in.) 


60 


ORIGINAL PAQf It 
OF POOR QUALITY 



Figure 5.6 Deflection versus time for two-layer plate ro*/90*] 
under uniformly distributed step load. 
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line In Mg. 5.8 Is for At ■ 0.1 x 10"^ sec, which gives slightly 
smaller deflections. 

5.3.3 Four-Layer Angle Ply (45V-45V45V-45* ) Cyl Indrical Shell 
Under Uniform Load 

Here we present results for a cylindrical shell which has the same 
geometry as problem 5.3.2. The four-layer angle ply (45*/-45*/45®/-45*) 
laminated shell Is singly supported on four edges and Is subjected to a 

A 

uniform step load P ■ 50. Fig. 5.9 contains a plot of the center 
deflection versus time for 2-0 and 3-0 elements. The two elements yield 
solutions that agree very well, except that the 2-0 element gives 
negative values of the deflection at the end of the cycle. The 
discrepancy Is due to the fact that the 2-0 element does not account for 
geometric changes from one time step to the next. 

5.3.4 Two-Layer Angle-Ply (45”/-45**) Spherical Shell Under Uniform 
Loading 

Consider a spherical shell with a ■ b ■ 10 In, R ■ 20 In and h * 

0.1 In, simply supported at four edges and excited by a uniform step 
load. The shell consists of two layers (45® /-45*). Figure 5.10 shows 

A A 

the center deflection versus time for P ■ 50 and P ■ 500 with time step 
0.2 X 10"^ sec. For the small load the curve Is re*'atively smooth 
conpared to that of the larger load. This Is due to the fact that the 

A 

geometric nonlinearity exhibited at P » 50 Is smaller conpared to that 

A 

at P * 500. 



I 

I Figure 5.8 Comparison of the center deflection obtained by 

j two different time steps for the problem discussed 

in Section 5.3.2. 

j 
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Figure 5.9 


Center xleflectlon versus time for four-1 ayei 
angle-ply [45*/-45®/45', -45*] cylindrical 
shell subjected to uniformly distributed loa 
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center deflection 
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— wxlO** for load parameter, P ■ 50 

% 

\ 

\ 

% 

— wxlO for load parameter, P ■ 500 

\ 

\ 

\ 

\ 

R/a ■ 2, a/b ■ 1, a/h ■ 100 


Time, txiO' 


Figure 5.10 Center deflection versus time for two-layer 
angle-ply [45*/ -45*] spherical shell under 
uniformly distributed step loading. 
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CHAPTER VI 

SUMMARY AND CONCLUSIONS 

6.1 Sunniary of the Present Study 

A special three-dimensional element based on the total Lagranglan 
description of the motion of a layered anisotropic conposit' medium Is 
developed, validated, and employed to analyze conposite shel Is. The 
element has the foliOMing options: 

• Geometrically linear and nonlinear analyses 

• Static and transient analyses 

• Natural vibration (linear) analyses 

• Plate and shell elements 

• Arbitrary loading and boundary conditions 

• Arbitrary lamination schemes and lamina properties 

The element can be used, with minor changes. In any existing general 
purpose program. 

6.2 Conclusions 

The present 3-0 degenerated element has conputatlonal simplicity 
over a fully three-dimensional element, such as those developed In [49- 
50], and the element accounts for full geometric nonlinearltles In 
contrast to 2-0 elements based on shell theories. As demonstrated via 
numerical exanples, the deflections obtained by the 2-0 shell element 
deviate from those obtained by the 3-0 element for deep shells. 

Further, the 3-0 element can be used to model general shells that are 
not necessarily doubly-curved. For exanple, the vibration of twisted 
plates cannot be studied using the 2-0 shell element discussed In 
Chapter 2. Of course, the 3-0 degenerated element is conputatlonal ly 
more demanding than the 2-0 shell theory element for a given problem. 
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In summary, the present 3-0 element Is an efficient element for the 
analysis of layered conposlte plates and shells undergoing large 
displacements and transient motion. 

6.3 Recommendations for Additional Study 

The 3-0 element presented herein can be modified to Include thermal 
stress analysis capability and material nonllnearltles. While the 
Inclusion of thermal stresses Is a sinple exercise, the Inclusion of 
nonlinear material effects Is a difficult task. An acceptable material 
model could be a generalization of the Ramberg -Osgood relation to a 
layered anisotropic medlun. Areas that require further study are the 
Inclusion of danping effects, which can be more significant than the 
shear deformation effects, and material damage effects. 

Acknowledgments 

The present study was conducted under a research grant from the 
Structures Research Section of NASA Lewis Research Center. The authors 
are very grateful for the support and encouragement by Or. C. C. Chamis 
of NASA/Lewis. It Is also a pleasure to acknowledge the typing of the 
manuscript by Mrs. Vanessa McCoy. 


68 


REFERENCES 


1. P. M. Naghdl, "A Survey of Recent Progrress In the Theory of 
Elastic 9ie11$.** Applied Mechanics Review Vol. 9. 1956, pp. 365. 

2. C. W, Bert, “Analysis of Shells," Structural Design and Analysis , 
Part I , C. C. Chamis (ed.), Vol. 7 of Composite Materials . L. J. 
Broutman and R. H. Krock (eds.). Academic Press, New York, 1974. 

3. S. A. Ambartsuiy an, “Calculation of Laminated Anisotropic Shells," 

Izvestia Akademiia Nauk Armenskoi SSR, Ser. Fiz. Mat. Est. Tekh. 
Nauk. , Vol. 6, No. '3ri9y3.'pp.T?: 

4. S. A. Ambartsinyan, Theory of Anisotropic Shells , Moscow, 1961; 
English Translation, NASA TT F-118, May 1964. 

5. S. B. Oong, K. S. Pister and R. L. Taylor, "On the Theory of 
Laminated Anisotropic Shells and Plates," Journal of Aerospace 
Sciences , Vol. 29, 1962, pp. 969-975. 

6. Y. Stavsky, “Bending and Stretching of Laminated Anisotropic 

Plates." Journal of the Enqi neerinq Mechanics Division, Proc. ASCE, 
Vol. 87, No. EM6, 1961, pp! 31. 

7. S. Cheng and B. P. C. Ho, “Stability of Heterogeneous Aelotropic 

Cylindrical Shells Uhder Combined Loading," American Institute of 
Aeronautics and Astronautics Journal, Vol. 1, No. 4, 1963, pp. 892- 
■ 35 ^: 

8. 0. E. Widera and S. W. Chung, “A Theory for Non-Homogeneous 
Anisotropic Cylindrical Shells," Journal of Applied Mathematics and 
Physics (ZA^f>), Vol. 21, 1970, pp. 378-399. 

9. S. T. Gulati and F. Essenburg, "Effects of Anisotropy in 
Axisymmetric Cyl indrical Shells," Journal of Appl led Mechanics, 

Vol. 34, Sept. 1967, pp. 659-666. 

10. J. A. Zukas and J. R. Vinson, "Laminated Transversely Isotropic 
Cylindrical Shells." Journal of Applied Mechanics. Vol. 38, 1971, 
pp. 400-407. 

11. J. M. Whitney and C. T. Sun, "A Refined Theory for Laminated 
Anisotropic, Cylindrical Shells," Journal of Applied Mechanics, 

Vol. 41, 1974, pp. 471-476. 

12. G. E. 0. Widera and 0. L. Logan, "Refined Theories for 

Nonhomogeneous Ani sotropic Cyl i ndrical Shells: Part I-Deri vation," 

Journal of the Engineering Mechanics Division Proc. ASCE, Vol. 106, 
No. Em 6, 7980, pp! 1053-137T: 

13. 6. E. 0. Widera and D. L. Logan, "Refined Theories for 

Nonhomogeneous Anisotropic Cylindrical Shells: Part II- 

Appl ication," Journal of the Engineering Mechanics Division Proc. 
ASCE. Vol. 106, No. EM6, 1980, pp. )07S-109b. 


69 


14. S. B. Dong, "Analysis of Laminated Shells of Revolution," Journal 
of the Engineering Mechanics Division Proc. ASCE, Vol. 92, EM6, 

1966, pp. 13S. 

15. S. B. Dong and L. G. Selna, “Natural Vibrations of Laminated 
Orthotropic Shells of Revolution," Journal of Conposlte Mateilals, 

Vol. 4, No. 1, 1970, pp. 2-19. 

16. E. A. Wilson and B. Parsons, "The Finite Element Analysis of 
Filament-Reinforced Axisymmetric Bodies," Fibre Science and 
Technology , Vol. 2, 1969, pp. 155-156. 

17. L. A. Schmit, Jr., and G. R. Monforton, "Finite Element Analysis of 
SandMich Plate and Laminated Shells with Laminated Forces," 

American Institute of Aeronautics and Astronautics Journal, Vol. 8. 

1970 , pp." 1454 - 1461 . 

t 

18. G. L. Thonpson and C. W. Bert, "Finite Element Analysis for Free 
Vibration of General Anisotropic Laminated Thin Shells," Conposite 
Materials In Engineering Design (Proc. 6th St. Louis Sym,'os1un, May 
ll-li, 1972), Norton, B. R. (ed.j, American Society for Metals 
Park, Ohio, 1973, pp. 72-80. 

19. S. C. Panda and R. Natarajan, "Finite Element Analysis of Laminated 
Shells of Revolution," Conputers and Structures, Vol. 6, 1976, pp. 

61-64. 

20. K. N. Shlvakunar and A. V. Krishna Murty, "A High Precision Ring 
Element for Vibrations of Laminated Shells," Journal of Sound and 
Vibration , Vol. 58, No. 3, 1978, pp. 311-318. 

21. K. P. Rao, "A Rectangular Anisotropic Shallow Thin Shell Finite 
Element," Conputer Methods In Applied Mechanics and Engineering, 

Vol. 15, 1978, pp. 13-83. 

22. P. Selde and P. H. H. Chang, "Finite Element Analysis of Laminated 
Plates and Shells," NASA CR-157106, 1978, 132 pp. 

23. A. Venkatesh and K. P. Rao, "A Doubly Curved Quadrilateral Finite 
Element for the Analysis of Laminated Anisotropic Thin Shell of 
Revolution," Conputers and Structures . Vol. 12, 1980, pp. 825-832. 

24. H. S. Hsu, J. N. Reddy and C. W. Bert, "Thermoelastlclty of 
Circular Cylindrical Shells Laminated of Bimodulus Conposite 
Materials," Journal of Thermal Stresses, Vol. 4, No. 2, 1981, pp. 

115-127. 

25. J. N. Reddy, “Bending of Laminated Anisotropic Shells by a Shear 
Deformable Finite Element, ‘‘ Fibre Science and Technology, Vol. 17, 
pp. 9-24, 1982. 


70 


26. G. Horrigmoe and P. G. Bergan, "Incremental Variational Principles 
and Finite Element Models for Nonlinear Problems," Computer ^thods 
m Applied Mechanics and Engineering . Vol. 7, 1976, pp, 201-217. 

27. W. Wunderlich, "Incremental Formulations for Geometrically 
Nonlinear Problems," Formlatlons and Algorithm In Finite Element 
Analysis, by K. J. Bathe, J. T. Oden and U. Wunderlich ed. , pp. 
193-239. 

28. A. Stricklin, W. E. Halsler and W. A. Von Riesemann, "Evaluation of 
Solution Procedures for Material and/or Geometrically Nonlinear 
Structural Analysis." AIAA Journal . Vol. 11, 1973, pp. 292-299. 

29. A. K. Ncor and S. J. Hartley, "Nonlinear Shell Analysis Via Mixed 
Isoparametric Elements," Computers and Structures . Vol. 7, 1977, 
pp. 615-626. 

30. T. Y. Chang and K. Sawami phakdl , "Large Deformation Analysis of 
Laminated Shells by Finite Element Msthod," Conputers and 

St ructures . Vol. 13, 1981, pp. 331-340. 

31. J. L. Sanders, Jr., "An Improved FI rst -Approximation Theory for 
Thin Shells," NASA Technical Report R-24, June 1959. 

32. C. H. Tsao, "Strains • Displacement Relations In Large Displacement 
Theory of Shells." AIAA Journal. Vol. 2. No. 11, 1964, pp. 2060- 
2062. 

33. J. N. Red<^ and W. C. Chao, "A Conparlson of Closed-Form and 
Finite-Element Solutions of Thick Laminated Anisotropic Rectangular 
Plates," Nuclear Engineering and Design, Vol. 61, 1981, pp. 153- 
167. 

34. J. N. Redct/, "On the Solutions to Forced Motions of Layered 
Composite Plates," AIAA Journal . Vol. 21, No. 4, 1983, pp. 

621-629. 

35. 0. C. Zienklewicz, R. L. Taylor and J. M. Too, "Reduced Integration 
Technique In General Analysis of Plates and Shells," International 
Journal for Numerical Methods In Engineering, Vol. 3, 1971, pp. 

m-m: 

36. B. Krakeland, "Nonlinear Analysis of Shells Using Degenerated 
Isoparametric Elements," Finite El erents in ^nllnear Mechanics. 
International Conference on Finite Elements In Nonlinear Solid and 
Structural Mechanics . Gello, Norway, Aug. 1977, pp. 265-284. 

37. K. J. Bathe, E. Ramm and E. L. Wilson, "Finite Element Formulations 
for Large Deformation Dynamic Analysis," International Journal for 
Numerical ^tethods In Engineering. Vol. 9 , 1975, pp. 353-386. 

38. G. A. Oupurls, H. D. Hibbit, S. F. McNamara and P. V. Marcal, 
"Nonlinear Material and Geometric Behavior of Shell Structures," 
Conputers and Structures. Vol. 1, 1971, pp. 223-239. 


71 


39. G. Horrigmoe and P. 6. Sergan, "Nonl i near Analysis of Free-Form 
Shells by Flat Finite Elements," tonputer Methods 1ji Applied 
Mechanics and Engineering . Vol. 16, 1978, pp. 1 1-35. 

40. S. Ahmad, B. H. Irons and 0. C. Zienklewicz, "Analysis of Thick and 
Thin Shell Structures by Curved Finite Elements," International 
Journal for Numerical Methods In Engineering, 2, 1970, pp. 419-451 . 

41. R. E. Nickel 1, "On the Stability of Approximation Operators In 
Problems of Structural Dynamics," Int. J. Solids and Structures, 
Vol. 7, pp. 301-319, 1971. 

42. R. E. N Ickell, "Direct Integration Methods In Structural 
Dynamics." J. Engng. Mech. Olv. ASCE . Vol. 99. No. EM2. 1972. 

43. G. L. Goudreau, and R. L. Taylor, "Evaluation of Numerical 
Integration Methods In Elastodynamics," J. Conputer Methods In 
Applied Mechanics and Engineering , Vol. 2, No. 1, pp. 69-97, 1973. 

44. T. B. Belytschko and B. J. Hsieh, "Nonlinear Transient Finite 
Element Analysis with Convected Coordinates," International Journal 
for Numerical Methods In Engineering. Vol. 7, pp. 255-27^, 1973. 

45. T. J. R. Hughes, and W. K. Liu, "Implicit-Explicit Finite Elements 

In Transient Analysis: Stability Theory, and Implementation and 

Numerical Examples," Vol. 45. No. 2. pp. 371-378, 1978. 

46. S. P. Chan, H. L. Cox and U. A. 8enf1eld, "Transient Analysis of 
Forced Vibrations of Conpl?x Structural Mechanical Systems," 

Journal of Royal Aeronant Soc. . Vol. 66, pp, 457-460, 1962, 

47. L. Fox and E. T. Goodwin, "Some New Methods for the Numerical 
Integration of Ordinary Differential Equations," Proc. Camb. Phil. 
Soc. Math. Phys. Scl. Vol. 45, 1949, pp. 373-388. 

48. 0. C. Zienklewicz, The Finite Element Method, 3rd ed. , McGraw-Hill 
Book Company. 

49. G. S, Ohatt, "Instability of Thin Shells by the Finite Element 
Method," IASS Symposium for Folded Plates and Prismatic Structures, 
Vienna, 1970. 

50. R. H. Leicester, "Large Elastic Deformations and Snap-Through of 
Shallow Doubly-Curved Shells," Ph.D. Thesis, Department of Civil 
Engineering, University of Illinois, Urbana, 1966. 

51. J, A, Stricklin, J. E. Martinez, J, R. Tlllerson, J. H. Hong, and 
W. E. Halsler, "Nonlinear Dynamic Analysis of Shells of Revolution 
by Matrix Displa -Tent Method," AIAA Journal. Vol. 9, No, 4, April 
1974, pp. 629-5 

52. S. Timoshenko and S. Wol nowsky-KrIeger, Theory o f Plates and 
Shells, McGraw-Hill, New York, 1959. 


72 


53. J. N. Red<t/. "Dynamic (Transient) Analysis of Layered Anisotropic 

Conposite-Materlal Plates," International Nwrical 

Methods In Engineering . Vol. 19, No. 2, pp. 237-25$, 

54. J. R. Dana and R. M. Barker, "Ihree-Olmenslonal Analysis for the 
Stress Distribution Near Circular Holes In Laminated Conposites," 
Research Report No. VPI-E-74.18, Virginia Polytechnic Institute and 
State University, Blacksburg, VA, 1974. 

55. J. N. Reddy and T. Kuppusany, "Analysis of Layered Conposite Plates 
by Three-Dimensional Elasticity Theory," Research Report No. VPI-E- 
82,31, Virginia Polytechnic Institute and State University, 
Blacksburg, VA, 1982. 


73 


ORIGINAL PAGE fs 
OF POOR QUALITY 


APPENODU. : element STIFFNESS ANO MASS MATRICES 

• AjjCS^^] ♦ (Ajg ♦ SBjg)(CS ^^3 ♦ ) 

^ (^66 * ^^ o ®66 ^ * ”7 

^1 

CK ^23 . (Ajg - SBjg)CS^h > (Agg - c20gg)CS^23T ^ 

+ (A^g + %B 2 g)[S^^] ♦ A^g [S °°3 

[K»] . (^ » ^ 12 )[S>»] * * c,(^ ♦ ^)] cs'O] 

- ^ * » 55 Cs“']) 

[k’^] • * (®16 * 

* (»66 * 'o» 56 )C 5 “] - CS'"’! 

Ck‘*] • ♦ * 16 ^ 5 ": * (Bj 5 ‘ 

* <»66 ‘ ^ A„ [S»»] 

[K^^] ■ AjjCS^^] ♦ (A35 - C^BjjKCS*^] ♦ Cs'^f ) 

* (*66 ■ ^ S ®66 * [ S °°] 


ck”: 


^®16 ■ Ao **16^^® ^ * ^*66 ■ ^ ^ 


* “JB CS''] - A„ [S»»] 

Cit”] • * (Bjj - C^DjjjCs'^] 


^ («66 ■ -_!ics“] 

«2 


74 


ORIGINAL PAGE W 
OF POOR QUALITY 




.22 


B2,CS“] . BjsCs'^] ♦ (8 


26‘ 


’26 


’ ^o°26 


12,T 


)CS“] 


♦ (8 


66 ’ ^o°66 


)CS^^] 


44 


[S°0] 




*45CS’^] 


AjgCS^^] + A 


>44 CS^*] • A 


45CS‘^]' 




* ^)CS»' 


] 



• 0„[s“] ♦ 0 j5([s‘ 2] * * AjjCs"] 

Ck'*] ■ 0|jCs‘2] ♦ DjgCs"! * . 0„[si2]^ 


* *45CS“] 

[k”] . OjjCS*^] . OjjCs"] * DjjCS^^] ♦ o^jCs'^]^ . ‘441:5°°] 
[M^^] = [M^^] = ^ Pj[S°°] 

[M^^] X [M^^] X P3[S°®] 

CM^^] = [M^^] -P- [S°°l P X P^Cl/R^rl/R^] 

^ , Vi'’sV’'l‘”'2 • (»•? ■ 0.'.2) > = si- . 0„ ■ I . 

Q a 


75 


ORfG^rjAL PAGE IS 
OF POOR QUALITY 


A 

«. P 
“i 

0 

'1 

»1J 


’IJkJl 


'ij 

lA} 


~T 

A 

e 

A# 

®1j 
{ F > 

♦i 

h 

[K] 

''i 

L 

[M] 


APPENDIX II; NOMENCLATURE 


- area of elenient 

- stretching stiffness matrix 

- dimensionless parameters of generalized acceleration 

- surface metric 

- strain in the i-th direction 

- strain in the ith direction at the reference plane 

- bend1ng>stretching stiffness matrix 

- configuration at time t 

- elasticity tensor 

- Eh^/12(1 - v^) 


- bending stiffness matrix 

- generalized displacement vector 
1 1 ^ 


- Young's modulus 

- unit vector along the global axes 

- unit vector along the local axes 

- linear incremental strain 

- curvilinear coordinate system 
•• force vector 

- rotations of the reference surface with respect tr 

- thickness 

- stiffness matrix 

- shear correction factor 

- nunber of layer 

- moment resultatnt 


- mass matrix 
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- Poisson's ratio 

- stress resultant 

- nonlinear incremental strain 

- inertias 

- distributed load 

- shear force resultant 

A A 

- rotations about unit vector ei.e-, 

- balance force vector 

- radi i of curvature 

- density 

- 2nd Piola Kirchhoff stress tensor 

- stress matri x 

- stress vector 

- length of line element of the shell 

- stress vector 

- interpolation function at node i 

- transformation matrix between displacement vector (u }and 
generalized displacement vector {A} 

- traction component 

- time increment 

- incremental displacement vector 

- displacement vector at time t 

- displacement gradient 

- vector along the local axes 

- volume 

- coordinate at time t 

- virtual work due to external loads 

- column of generalized nodal displacements 
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